1 



oo 

On 
On 



O 

O 
o 

(N 



> 
O 

o 
o 

00 
On 

D 



X 



FERMILAB-Conf-98/333-T 
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We estimate statistical and systematic errors in the extraction of the CKM ratio |14d|/| V cs | from exclusive 
D-meson semileptonic decays using lattice QCD and anticipated new experimental results. 



1. INTRODUCTION 

High statistics experimental studies of D 
mesons, such as E831 (FOCUS) and E791, are ex- 
pected to yield O(10 6 ) fully reconstructed D de- 
cays with better momentum resolution than pre- 
ceding experiments [1,2]. Anticipating high pre- 
cision experimental results, we consider how well 
the CKM ratio | V c d | / 1 V cs | can be determined from 
a ratio of D — > nev to D — > Kev decay rates and 
lattice QCD. 

For the Cabibbo allowed D — > Kev decay, 
lower statistics and poorer resolutions in ear- 
lier experiments meant hadronic matrix elements 
could be described by functions of only two pa- 
rameters: the CKM matrix element times a form 
factor at zero momentum transfer and the pole 
mass in a pole dominance Ansatz. The PDG 
quotes a 3% error for |Ks|/+^(0) based upon such 
techniques [3]. The uncertainty due to this Ansatz 
is difficult to quantify. They quote a 14% uncer- 
tainty in the theoretical determination of f? (0) 
resulting in a 15% uncertainty for \V CS \. 

\V c d\ is determined with an error of 7% from v- 
v production of charm off of valence d quarks [3]. 
Combining this uncertainty and that for \V CS \ we 
deduce an error of about 17% in the current best 
value of \V cd \/\V cs \- 

We examine uncertainties in determining 
I Vcd | / 1 Vcs | directly from a ratio of semileptonic 
decay rates. We expect that both experimen- 
tal and theoretical systematic errors common to 
D — > Kev and D — > itev will be reduced in ratios. 
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2. PROCEDURE 

We suggest the ratio of CKM matrix elements 
can be determined by a double ratio 

I cut \ 
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The maximum recoil momentum cut, |p| c "*, will 
be adjusted to minimize combined experimental 
and theoretical errors [4]. The ratio of partial 
widths 
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is to be determined experimentally. T^/ K (\p\ cut ) 
is provided by theory calculation. In the D-meson 
rest frame, 
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where hadronic matrix elements are expressed 
conveniently by form factors f¥. 
Matrix elements 

(X(p x )\V^\D( PD )) 

y/2Ejj2Ex 

are determined in the continuum limit of our lat- 
tice results. do not depend upon the heavy 
quark mass in the infinite mass limit. The con- 
tinuum relation 
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(PD - P x)J X (q 2 )] /V2E D 2E X (5) 
is used to extract the form factors. 
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Figure 1. Decay rates (with statistical errors) ver- 
sus daughter recoil momentum. 

3. MATRIX ELEMENTS 

We use the tadpole-improved Shcikholcslami- 
Wohlert quark action with the plaquette- 
determincd mean field coefficient uq. Wc obtain 
O(a)-improved matrix elements using the Fermi- 
lab interpretation for heavy quarks and tree-level 
O(a)-improved currents [5]. Note that the vec- 
tor current renormalization factor cancels in ratio 

Ttt/K- 

We have computed three-point correlators at 
(3 = 5.7, 5.9 and 6.1. We determine the lattice 
spacing, a, using the charmonium IP-IS splitting 
and the bare charm mass using the spin-averaged 
15 charmonium kinetic mass. The same gauge- 
field ensembles and procedures were used to study 
B and D meson decay constants [6]. Reference [7] 
describes our B — > nev results on these lattices 
and provides more details on our three-point func- 
tion calculations. 

For illustration, and to reduce some known sys- 
tematic errors, we take |p| cu * = 0.7 GeV. For this 
preliminary study we take our best estimate for 
matrix elements from (3 — 5.7. We approximate 
the D — > Kev decay by D s — ► ~K s -^ev where a 
pseudoscalar "7r S s" meson consists of s and s va- 
lence quarks. We denote the ratio in Eq. 3 for the 
D s decay by T n ^ S g. The analysis for D — > Kev is 
underway 

In Fig. 1 wc show differential decay rates for 
(3 = 5.7. The ratio of areas for |p| cut = 0.7 GeV 
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Figure 2. Comparison of MJf(p S s, 0) for the de- 
cay of B s and D s mesons at (3 = 5.7. 



gives T-x/ss = 1.53 ± 0.27 where the error is sta- 
tistical. Systematic errors are the subject of the 
rest of the paper. 

We use all three lattices to study the cutoff de- 
pendence of D s — > TTsjeiy matrix elements. We 
compare partial widths at (3 = 5.7 to partial 
widths computed from M. 8 ^ extrapolated to the 
continuum linearly with the lattice spacing. Ma- 
trix elements A4™(p,0) for a given p on each 
lattice were found by interpolating among M.'jf 
values calculated for discrete values of p. Dis- 
cretization errors in the width increase with p, as 
expected, reaching 21% for |p| = 0.7 GeV. Mo- 
mentum dependent errors of this magnitude were 
anticipated by errors estimates for free-field tree- 
level quark matrix elements [4]. Errors for the 
pion width are expected to be similar to those 
for the ss width. These errors, however, will tend 
to cancel in ratios. Hence we estimate residual 
cutoff errors are below 10% for T ff / S j(0.7). 

Our bare charm masses were determined in 
quarkonia. Charm masses determined using 
heavy-light mesons differ due to discretization er- 
rors and the quenched approximation. We find 
at most a 10% difference in bare m c compar- 
ing the two methods. Figure 2, which com- 
pares _M* s (p, 0) calculated for D s and B s mesons, 
shows that I/tuq corrections are small in this 
case [4]. Numerical results and pole dominance 
arguments indicate that M. 1 * has somewhat more 
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1/rriQ dependence than M s ^ s . We estimate the 
uncertainty in adjusting m c leads to a 1% error 
in T w/s3 . 

The quenched approximation leads to system- 
atic differences in lattice spacing determinations. 
We take the charmonium IP-IS lattice spacing as 
our best estimate of a. Repeating our analysis for 
another reasonable choice for the lattice spacing, 
a(/ftr), determined by the kaon decay constant, 
changes T^/ S g by about 1%. This almost certainly 
underestimates the error due to quenching. 

Correlators are computed for light quark 
masses between about 0.4 to 1.2 times m s . For 
lighter masses, correlators begin to show in- 
creasing fluctuations caused by nearby unphysi- 
cal Dirac operator zero-modes found on a small 
subset of our gauge configurations. Hence we 
must obtain matrix elements for up and down 
quarks by "crural" extrapolations of matrix ele- 
ments computed for masses close to the strange 
quark. Near zero recoil momentum these extrapo- 
lations are more difficult because matrix elements 
are sensitive to the light quark mass. This can 
be understood as the effect of nearby poles in a 
pole dominance picture. The effect upon T^i S g 
is small, however, since the decay rate is kine- 
matically suppressed by p 4 (see Eq. 3). Matrix 
elements show less dependence upon quark mass 
at larger recoil momenta. Statistical noise in- 
creases with recoil momenta, however, making 
extrapolations there noisier. Hence we find a 
"window" at moderate recoil momentum where 
chiral extrapolations are most dependable. By 
varying our extrapolation procedures we estimate 
chiral extrapolations contribute a 10% error to 
T„/ S g(0.7). If correlators could be reliably com- 
puted for masses significantly below m s our chi- 
ral extrapolation uncertainties would decrease. 
The Modified Quenched Approximation shows 
promise in this respect [8]. 

Errors due to varying fitting procedures and 
residual excited state contamination are esti- 
mated to be about 2%. 



4. CONCLUSIONS 

Table 1 



Table 1 

Summary of errors in \V c d / \V CS \ at (3 = 5.7 


source 


% 


M.C. statistics 


18 


fits, excited state contamination 


2 


residual cutoff dependence 


< 10 


a determination 


1 


bare m s det., extrapolation to mi 


10 


bare m c determination 


1 



nary results are based upon our (3 — 5.7 analy- 
sis. We are completing the analysis of K and ir 
decays on our finer [3 = 5.9 lattice and continu- 
ing to study cutoff effects using our 6.1 lattice. 
Smaller a will result in smaller cutoff effects than 
our (3 = 5.7 results but with similar statistical er- 
rors. The error due to the quenched approxima- 
tion has yet to be fully investigated. A quenching 
error of 10-20% would not be unexpected, how- 
ever, in the deviation of T^/k from unity. 

Theoretical uncertainties in |V c d|/|V cs | have 
been estimated for our procedure with |p| cu * = 
0.7 GeV. We selected this cut to eliminate the 
largest momentum dependent errors and chiral 
extrapolation errors. Before quoting a value for 
|14d|/|Ks|, |p| cu * should be chosen to minimize 
the total theoretical and experimental errors. 
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